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The natural thermoluminescence (TL) of meteorites reflects 
their irradiation and thermal histories. Virtually all ordinary 
chondrites have been irradiated long enough to have reach 
saturation natural TL levels, and thus natural TL levels in 
these meteorites are determined largely by thermal history. The 
primary heat source for most meteorites is the Sun, and thus 
natural TL levels are determined primarily by the closest ap- 
proach to the Sun, i.e., perihelion. By converting natural TL 
levels to perihelia, using an assumed albedo typical of meteoroid 
bodies, it is found that most ordinary chondrites had perihelia 
of 0.85 to 1.0 AU prior to reaching Earth. This range is similar 
to that calculated from meteor and fireball observations. All 
common classes of ordinary chondrites exhibit similar perihelia 
distributions; however, H and LL chondrites that fell in the 
local morning differ in their natural TL distribution from those 
that fell in the local afternoon and evening. This is consistent 
with earlier suggestions that time of fall reflects orbital distribu- 
tion. The data also suggest that the orbits of some of the H 
chondrites cluster and may have come from a debris “stream” 
of meteoroids. If meteorites can exist in “orbital groups,” sig- 
nificant changes in the types and number of meteorites reaching 
Earth could occur on the < 10 5 -year time scale. 1997 Academic Press 


INTRODUCTION 

Most meteorites are thought to be samples of asteroids 
(Wetherill 1985); however, it is unlikely that meteorites 
are random representative samples of their parent bodies. 
Their passage to Earth depends on their original orbit and 
the stochastic nature of the collisional events that ejected 
them into space (Wetherill 1985, Greenberg and Nolan 
1989, Binzel el al. 1993). It has been proposed that some 
fraction of the meteorite flux derives from coherent orbital 
“streams" of meteoroids (Hallidav 1987, Dodd etal. 1993). 
This idea is indirectly supported by orbital clustering of 
meteor-producing bodies and temporal clusters of meteor- 
oid impacts on the Moon (Olsson-Steel 1988, Halliday et 
al. 1990, Drummond 1991, Oberst and Nakamura 1991). 
There is. however, little evidence that, most meteorites 
were members of meteoroid streams at the time of their 
fall (Wetherill 1985). 


The possible existence of meteoroid streams, even as 
relatively minor contributors to the terrestrial flux of extra- 
terrestrial material, is important in understanding the or- 
bital evolution of meteoroid bodies, the rarity of potential 
asteroidal sources for ordinary chondrites, which account 
for more than 60% of meteorite falls (Greenberg and Nolan 
1989), and apparent differences between meteorites cur- 
rently failing and those that fell over the last few hundred 
thousand years and have been stored in the ice sheet of 
Antarctica (Koeberl and Cassidy 1991, Benoit and Sears 
1993a, Wolf and Lipschutz 1995). 

Precise orbits are known for only 4 meteorites (Wetherill 
and Chapman 1988, Brown et al. 1994). Orbital parameters 
have been estimated for another 40 meteorites from eye- 
witness accounts of their fall (Simonenko 1975, Wetherill 
and Chapman 1988), but these are of low accuracy. Limited 
orbital information can also be obtained from the local 
time of fall (Wetherill 1968). Meteorites with perihelia -1 
astronomical unit (AU) should usually strike the trailing 
side of the Earth, or the local pm, while meteorites with 
perihelia <1 AU should be more evenly spread over the 
leading and trailing sides. The large abundance of ordinary 
chondrites falling in the local pm relative to am is interpre- 
ted as indicating that most of these meteorites had perihelia 
of -1 AU. 

In this study we use the natural thermoluminescence of 
individual meteorites to estimate their closest approach to 
the Sun. Natural thermoluminescence is energy stored in 
crystals of certain minerals by ionizing radiation, such as 
high-energy galactic cosmic rays (Aitken 1985, McKeever 
1985). This energy can be released in the form of visible 
light by heating. 

The natural thermoluminescence (TL) of meteorites re- 
flects both their thermal and irradiation history (McKeever 
1980. Benoit et ciL 1991). Virtually all ordinary chondrites 
have been exposed to galactic cosmic radiation for millions 
of years (Crabb and Schultz 1981), far longer than the 
~10 5 years required to obtain TL equilibrium (Benoit et 
al. 1991). Depth effects on TL levels are minimal, TL 
profiles generally exhibiting <30% variation in even the 
largest meteorites (Benoit and Chen 1996). The intense 


281 

0019*1035/97 S25.00 
Copyrighl £ 1997 by Academic Press 
All rights of reproduction in any form reserved. 



282 


BENOIT AND SEARS 


heating of atmospheric passage drains natural TL in the 
outermost few millimeters of meteorites (Sears 1975, Be- 
noit and Chen 1996). but we have deliberately avoided 
sampling these regions in this study. Therefore, the natural 
TL of modern falls is controlled largely by the maximum 
temperature experienced in space, which occurs at perihe- 
lion. The only exceptions are very large meteoroid bodies 
and meteoroids that have experienced shock heating within 
the last HP years, and both these exceptions appear to be 
very rare among the ordinary chondrites. 

Following the treatment of Garlick and Gibson (1948). 
natural TL can be considered as the filling of “traps" by- 
electrons put into the conduction band as a result of inter- 
action with energetic particles, such as high-energy cosmic 
rays. The exact nature of these “traps" is presently un- 
known. but is generally considered to involve crystallo- 
graphic defects and impurities (Prescott and Fox 1993). In 
crystals, populations of these “traps" are the source of TL, 
the electrons being freed from the traps and emitting a 
photon of light as the sample is heated. A typical TL 
experiment resulting in a “glow curve" involves heating a 
sample up to about 500°C in an inert atmosphere at a 
constant rate, while monitoring the amount of light given 
off by the sample as a function of temperature. In the 
present study we rearrange the basic equilibrium TL equa- 
tion to solve for temperature (McKeever 1980, Benoit et 
ai 1991), 


T = 



-E 

(\N/n] - 1 ) ln(2)r 
sR 


( 1 ) 


where n/N is the measured TL level relative to saturation. 
R is the radiation dose required to fill half the TL “traps," 
k is Boltzmann's constant, and r is the radiation dose rate. 
We use a dose rate of 0.01 Gv/year (1 rad/year) in our 
calculations; the actual dose rate experienced by samples 
will depend on the depths of the samples in the meteoroid, 
a factor w'e discuss below. E and ,v are characteristic of 
the TL phosphors, sodic feldspar in the case of ordinary 
chondrites, at a particular portion of the glow curve and 
they are essentially constants for ordinary chondrites 
(McKeever 1980). We can then calculate the distance from 
the Sun (X) from the calculated temperature assuming 
blackbodv behavior. 


0.25a// 
e(tT 4 


( 2 ) 


where a and e are the absorptivity and emissivity of the 
surface of the meteoroid body, El is the solar constant, and 
a is the Stefan-Boltzmann constant (Melcher 1981). Using 



FIG. 1. Equilibrium natural thermoiuminescence of ordinary chon- 
drites at 250°C in the glow curve as a function of closest approach to the 
Sun (perihelion). Perihelion is given in astronomic units (AU). 


a blackbodv approximation in this case is appropriate be- 
cause no asteroids, including those observed closely by 
space probes (Helfenstein et ai 1994). have atmospheres 
and the abundances of radioactive elements are too low in 
ordinary chondrites for these bodies to have any significant 
degree of internal heating during the last few million years. 
As a simplifying assumption, we set e to unity and estimate 
a value of a of 0.85 from the albedos of meteorite samples 
in the laboratory and remote-sensing measurements of as- 
teroid surfaces (Chapman and Salisbury 1973. Helfenstein 
et ai 1994). In Fig. 1 we show the relationship between 
natural TL level and perihelia for ordinary chondrites de- 
rived from Eqs. (1) and (2). This calculation shows that 
equilibrium natural TL levels exhibit a strong dependence 
on perihelia for perihelia <1.1 AU. We show curves for 
three dose rates, ranging from 0.005 to 0.05 Gy/year. In 
the following discussion we quote perihelion values using 
the 0.01 -Gv/year curve, an apparently typical value for 
ordinary chondrite meteoroid bodies (Benoit and Chen 
1996. and references therein). The 0.005- and 0.05-Gy/year 
curves are more appropriate for samples from very small 
and very large meteoroid bodies, respectively, and it is 
apparent that differences in meteoroid size have fairly min- 
imal effect on our interpretation. 

When we give a perihelion value for a particular meteor- 
ite, numerous uncertainties are included, some of which 
cannot be estimated with accuracy. Equation (2) is suscep- 
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tible to uncertainties in albedo. Meteorite samples do ex- 
hibit some range of albedoes (Chapman and Salisbury 
1973) and the relationship between the albedoes of meteor- 
ites and those of the surface of meteoroid bodies, of prime 
importance in this study* is not known. The major TL 
parameters. E and 5 , appear to exhibit little or no significant 
variation among ordinary chondrites (e.g., McKeever 
1980), but, as mentioned above, some uncertainty is im- 
posed because the size of the meteoroid body before atmo- 
spheric passage and the position of the TL sample in the 
body are generally not known, resulting in some uncer- 
tainty in dose rate. Because these uncertainties cannot be 
quantified, we base most of our analysis on our original 
TL data. The quotation of perihelion value, either below 
or in our figures, should be regarded as an interpretation 
of these data. 


EXPERIMENTAL 

For each meteorite we crushed an approximately 250- 
mg chip to a fine powder and removed metal (which is 
nonluminescent) with a hand magnet. We placed 4-mg 
aliquots of powder in copper pans which were then heated 
in a nitrogen atmosphere at a rate of 7.5°C/sec from room 
temperature to 500°C in a Daybreak Nuclear and Medical 
Systems TL apparatus fitted with blue bandpass and IR 
filters (Corning 7-59 and 4-69). Three aliquots of each 
sample were analyzed and the results averaged. Induced 
TL was measured by the same procedure 5 min after expo- 
sure to a 20-Gy dose from a ^Sr-^Y source. We used the 
ratio of the low-temperature (250°C in the glow curve) 
peak to the high-temperature (400°C) peak to find the 
equivalent dose at 250°C in the glow curve, using a pre- 
viously determined laboratory conversion equation (Hasan 
et al. 1987). The equivalent dose is the amount of radiation 
exposure needed to duplicate the observed natural TL 
level at a given temperature in the TL glow curve. A sample 
of a laboratory standard (the Dhajala H3.8 meteorite) was 
analyzed at the start and end of each measurement period 
as both a normalization standard and a check on the long- 
term stability of the apparatus. 

Uncertainties in the natural TL data from analytical 
sources, as determined by triplicate measurements on each 
sample, are quite small, being generally a few percent of 
the natural TL level, except for TL levels below about 15 
Gy that may have analytical uncertainties of up to 10-15%. 
We did not collect multiple samples from each meteorite 
and therefore we cannot directly assess the uncertainties 
caused by depth dependence; however, detailed TL studies 
of other meteorites indicate that typical sample-to-sample 
heterogeneity in meteorites introduces uncertainties of no 
more than 10% (Benoit and Chen 1996). With the excep- 
tion of meteorites with very low natural TL levels ( < 15 


16 

14 

12 

10 

8 

6 

4 

2 


(0 

0) 

o 

c 

£ 

L_ 

3 

o 

u 

o 


-O 

E 

3 


h L Chondrite 


0.6 
1 — 


Perihelion (AU) 

0 7 0 8 0.9 

1 I I — 


1.0 


H Chondrite 


J. 


■ ■■ 




II. 


n 


T 1 “ 1 

LL Chondrite 


n M n, 


CL 


J JS 

h H Chondrite "Stream 1 B" 


J.1 


1 10 100 1000 
Natural Thermoluminescence (gray) 


FIG. 2. Natural thermoluminescence at 250°C in the glow curve for 
ordinary chondrites. The distributions for H, L, and LL chondrites are 
fairly similar and correspond to orbits with perihelia between -0.75 and 
- 1.0 AU. Members of a possible H chondrite orbital debris stream ( Dodd 
et al, 1993) have a distribution similar to that of other ordinary chondrites. 


Gy), uncertainties in our data are no larger than the size 
of width of the boxes on our histograms (Fig. 2). 

RESULTS 

We have determined natural TL levels for 118 equili- 
brated ordinary chondrites (Fig. 2). Our previous work 
indicates that the natural TL levels of equilibrated ordinary 
chondrites are not affected by nonthermal losses or anoma- 
lous fading (Wintle 1973, Sears et al. 1991). There are no 
significant differences in natural TL between the major 
classes, the H. L, and LL chondrites, all classes having 
natural TL values of 50 to 1000 Gy, with a few outliers at 
lower values. 


DISCUSSION 

We can convert the natural TL data for individual mete- 
orites to perihelion using Eqs. ( 1 ) and (2). As noted above, 
such a conversion carries uncertainties that cannot be 
quantified, but differences of a few hundredths of an AU 
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FIG. 3. An interpretation of the natural TL data for ordinary chon- 
drites ( Fig. 2), using Eqs. ( 1 ) and (2) to convert the TL data to perihelia. 
For comparison, perihelion distributions derived from chondrite-like 
fireballs ( Wetherill and ReVelle 1981 ). visual radiants of observed mete- 
orite falls (Simonenko 1975). and MORP fireballs ( Hallidav ei ai 1983) 
are also shown. In general, the distribution derived from TL is similar 
to other distributions, and is most similar to those exhibited by chondrite- 
like fireballs and the distribution derived from MORP fireballs. 


are probably significant in comparisons between meteor- 
ites. The range of natural TL levels corresponds to perihe- 
lia of 0.75 to 1.0 AU, with an average around 0.95 AU, 
with just a few outliers with perihelia <0.75 and >1.0 AU 
(Fig. 2). This distribution is in general agreement with the 
distribution suggested by observed meteorite falls and also 
with the more numerous meteor and fireball radiant data 
(Fig. 3) (Wetherill and Chapman 1988). The distribution 
derived from the natural TL of meteorites is most similar 
to the distribution observed for fireballs attributed to ordi- 
nary chondrite objects, defined by two different methods 
(Wetherill and ReVelle 1981, Hallidav et ai . 1983), and 
differs from the distribution derived from visual radiants 
of ordinary chondrite falls by having a larger proportion 
of meteorites with perihelia between 0.8 and 0.9 AU. The 
most likely cause of this difference, as discussed by Wether- 
ill and ReVelle (1981), is biases inherent in the estimates 
of orbits from radiants observed by untrained persons. 


The outliers with natural TL levels <50 Gy (Fig. 2) are 
not discussed further here and are not considered in our 
statistics, because it is possible that some of these meteor- 
ites had their TL levels lowered by recent shock heating 
or we may have accidently sampled material near the fusion 
crust created during atmospheric passage. Thus, for these 
meteorites, the natural TL level could reflect a heating 
event other than closest approach to the Sun. 

As was noted above, one of the key observations that 
led to the delineation of the orbits and source of meteorites 
was the am/pm time-of-fall distribution (Wetherill 1968). 
We would expect such differences also to be reflected in 
the TL data and we do find such differences in at least two 
of the major ordinary chondrite classes (Fig. 4). The H 
chondrites that fell in the local am have a mean natural 
TL level similar to that of chondrites that fell in the pm 
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FIG. 4. Natural thermoluminescence distributions for H. L. and LL 
chondrites that fell in the am and pm. local time. For H chondrites the 
am falls exhibit a lower degree of variation than pm falls, whereas for 
LL chondrites the am falls generally have lower TL levels than pm falls. 
Onlv meteorites with natural TL levels >50 Gy are considered in our 
statistics and the total number used for each class and time-of-fall is given 
by n. 
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(about 300 Gy), but the am falls had a much smaller stan- 
dard deviation, lacking the large number of samples with 
natural TL of 100-300 Gy present in the pm falls (Fig. 2). 
The F test, used to examine equality of variance between 
data sets (Davis 1986), indicates that this difference is sig- 
nificant at the 95% confidence level and is thus probably 
not simply due to the smaller number of am falls in our 
study. The LL chondrites that fell in the local am have 
lower natural TL levels than PM falls, the difference being 
significant at the 99% confidence level according to the t 
test, but the distributions have similar variances. There 
are no significant differences in the TL distributions of L 
chondrite am and pm falls. 

Our interpretation of these data is that the am falls of 
LL chondrites have perihelia between 0.85 and 0.90 AU. 
whereas the pm falls are dominated by meteorites with 
perihelia between 0.85 and 1.0 AU. On the basis of the 
orbital dynamics arguments of Wetherill (1968), it can be 
argued that these data support the conclusion that LL 
chondrites are a steady-state population derived from the 
asteroid belt via orbital resonances that preferentially place 
them in orbits with perihelia of —1 AU. These bodies are 
short-lived, being removed from the inner solar system by 
interactions with planets on the 10 6 -year time scale, and 
only a few objects achieve orbits with perihelia significantly 
<1 AU. Our data indicate that very few LL chondrites 
that reach Earth have perihelia <0.8 AU. An alternative 
interpretation of these data is that the am and pm fall 
samples experienced different dose rates as a result of am 
falls being larger than pm falls. Our investigations of TL 
profiles in meteorites, however, suggest that the observed 
difference would require all am falls to be larger than 
typical meteoroid bodies (Benoit and Chen 1996). Such a 
conclusion is not supported by the size-dependent ratios 
of cosmogenic noble gases in LL chondrites (Graf and 
Marti 1994). 

The H chondrites have a more complicated apparent 
orbital distribution (Fig. 4). While the am falls have perihe- 
lia values of —0.95 AU, the pm falls, which should be 
dominated by meteorites with perihelia of about 1 AU, 
exhibit a range of 0.80 to 1.0 AU with a significant number 
<0.9 AU. One possible interpretation of this result is that 
the pm falls exhibit an extreme range of preatmospheric 
sizes and the am falls are very homogeneous in size. This 
interpretation is not supported by the ratios of cosmogenic 
noble gases found in these meteorites such as 2, Ne and 
“Ne, which vary as a function of meteoroid size (Eugster 
1988, Michel et ai 1991). 

An alternative interpretation of these data is that the H 
chondrites with TL levels between 100 and 300 Gy in the 
pm falls experienced shock heating within the last 10 5 years, 
the time required for TL to rebuild to thermally controlled 
equibrium: however, these meteorites display no petro- 
graphic evidence for high levels of shock. Another possibil- 


TABLE I 

H Chondrites in an Inferred Orbital Cluster 


Meteorite 

Class 

Cosmic ray- 
expos u re 
age ( Ma) 

:: Ne/ :i Ne 

Barbotan 

H5 

9 

1.10 

Beaver Creek 

H4 

29 

1.04 

Charsonville 

H6 

55 

1.06 

Forest Citv 

H5 

6 

1.08 

Gopaipur 

H5 

10 

1.24 

Lillaverke 

H5 

— 

— 

Noblesville 

H4 

— 

— 

Pokhra 

H5 

3 

1.24 

Pribram 

H5 

15 

1.11 

Rochester 

H6 

25 

1.07 

St. Germain du Pinel 

H6 

36 

1.14 

Zhovtnevy 

H5 

42 

1.08 


ity is that these meteorites were in orbits with perihelia of 
about 0.8 Au but moved to orbits with perihelia of — 1 AU 
within the last 10^ years. Such rapid orbit changes are 
possible, or even expected, for Earth-crossing asteroids as 
a result of orbital interactions with planets (Galibina and 
Terent’eva 1987, Farinella etal 1993). We have previously 
described a large Antarctic meteorite shower that appar- 
ently underwent rapid evolution from orbits with perihelia 
>1.2 AU to —1 AU (Benoit and Sears 1993b). This inter- 
pretation, however, requires that these meteorites changed 
orbit as a group. This could occur if these meteorites were 
part a single body that fragmented after changing orbit, 
but this seems unlikely due to their petrographic diversity 
(encompassing petrographic types 4, 5, and 6) and wide 
range of cosmic ray exposure ages (Table I). Alternatively, 
they could have formed an orbital debris ‘‘stream,” similar 
to those that produce meteor showers (Olsson-Steel 1988. 
Halliday et at. 1990). In this case, these meteorites reach 
Earth only when Earth's orbit intersects the stream. 

Dodd et at (1993) have examined the time-of-fall statis- 
tics for H chondrites in an attempt to identify possible 
orbital streams among historical falls. They proposed that 
a small subset of H chondrites (including some of the pm 
falls with apparent perihelia of —0.8 AU in our data set) 
were in an orbital stream, called H Cluster 1. Our analysis 
of 7 of 11 potential members of this stream, however, 
indicates that there is no difference between their TL distri- 
bution and that of other H chondrites (Fig. 2), suggesting 
that these meteorites did not have similar perihelia in their 
recent history and thus were probably not members of a 
modern orbital stream. This conclusion, however, rests on 
the assumption that a meteoroid stream has a stable orbit 
on the time scale reflected in the collection of modern falls, 
namely, on the scale of a few decades up to about a century, 
and thus that such streams are similar to the debris streams 



286 


BENOIT AND SEARS 


50 


-o 

to 

L_ 

8 '° 

§ 5 

(/) 

0) 

c 

E 

3 

O 

E i 

L— 

Q) 

-C 

t 0 5 
TO 
w 
_3 

TO- 

z 

0 1 


FIG. 5. Natural thermoluminescence levels of H chondrites from the 
proposed Cluster IB orbital stream of Dodd et at. (1993). plotted as a 
function of their year of fall. Meteorites in triangles along the bottom 
axis were not analyzed in this study. Meteorites that fell before about 
the year 1870 exhibit a fairly tight range of natural TL levels, consistent 
with perihelia slightly less than 0.95 AU. Those that fell after the year 
1885 have lower natural TL levels, possibly reflecting perihelia less than 
about 0.85 AU. AmN. Ambapur Nagla; BeC. Beaver Creek; CrR, Cross 
Roads; Gop, Gopalpur; Ker, Kernouve; Oto. Otomi: Pok. Pokra; Sea, 
Searsmont; Sla. Slavetic; Tyl, Tysnes Island. 
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responsible for meteor showers. Dodd et ai (1993) point 
out that a subset of the proposed H Chondrite Cluster 1, 
the IB stream, is not so well behaved and, although the 
members of the proposed stream fell over a period of about 
30 years (1865-1895), there is a pronounced gap within this 
range of about 12 years (1872-1883) in which no potential 
members fell. They suggest this reflects the narrow width 
of the stream and perhaps an orbit inclined to the ecliptic. 
We note that there is an apparent difference in the TL 
levels of meteorites from the proposed stream depending 
on their year of fall (Fig. 5). Meteorites that fell before 
the gap appear to have moderate natural TL levels, consis- 
tent with perihelia of slightly less than 0.95 AU. Meteorites 
that fell after the gap, however, have significantly lower 
natural TL levels, consistent with perihelia of <0.85 AU, 
and one, Assisi, may have had a perihelion of about 0.6 
AU (Benoit et ai 1991). It is unlikely that Assisi experi- 
enced TL loss due to a recent shock event, having a low 
petrographic shock level (Dodd et ai 1993) and a TL sensi- 
tivity of 4,3 ± 0.2 relative to the Dhajala meteorite, similar 
to other unshocked equilibrated ordinary chondrites (Haq 
et ai 1988). The apparent gap in the proposed Cluster IB 
could reflect a change in the orbit of the stream, with the 
members of the stream achieving much lower perihelia 
after the gap, or. conversely, this could reflect the sampling 


of different portions of a narrow stream by slight variations 
in the Earth's orbit. From the perspective of Fig. 4, in 
contradiction of the summary data of Fig. 2, it appears 
more likely that the proposed Cluster IB is actually a 
meteoroid stream, albeit a rapidly evolving one. with the 
stream breaking up in the second encounter. We stress, 
however, that this conclusion is of necessity based on a 
very small number of meteorites. 

CONCLUSIONS 

Only a small fraction (<25%) of the FI chondrites we 
have analyzed are involved in an inferred orbital cluster. 
In general, our data indicate that the remainder fit into 
the model of meteorid orbital distribution suggested by 
Wetherill (1968). Most ordinary chondrites appear to have 
been in orbits with perihelia between 0.8 and 1.0 AU prior 
to reaching Earth, and therefore they were probably de- 
rived as individuals from the asteroid belt through orbital 
resonance with the major planets. We find differences be- 
tween the orbital distributions of meteorites falling in the 
local morning and those falling in the afternoon and eve- 
ning. As predicted from indirect data, most ordinary chon- 
drites were in orbits with perihelia >0.85 AU, but we 
suggest that a small fraction of meteorites have experi- 
enced recent histories involving rapid orbital changes. We 
find some support for the existence of a rapidly evolving 
stream of FI chondrites, previously suggested on the basis 
of time-of-fall data. In general, it can be argued that minor, 
but perhaps important, changes have occurred in the num- 
bers and types of meteorites reaching Earth on the <10 5 - 
vear time scale. 
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